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Abstract 

The time evolution of a class of completely integrable discrete Lotka-Volterra system is 
shown not unique but have two different ways chosen randomly at every step of generation. 
This uncertainty is consistent with the existence of constants of motion and disappears 
in both continuous time and ultra discrete limits. 

1 Introduction 

Generally speaking we are not able to specify all of solutions of a given equation. A 
deterministic equation is either completely integrable or nonintegrable. In case of a non- 
integrable system, some of orbits converge into stable orbits, otherwise behave chaotic. 
Moreover, in spite of the deterministic nature of equation, chaotic orbits are very sensitive 
on their initial conditions. Namely one can not predict behaviour of orbits which had 
infinitesimally small difference of initial values. 
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Soliton equations are, therefore, distinct among other nonlinear equations. They are 
completely integrable under wide class of initial data. Solutions are often given explicitly 
in terms of r functions. One of important features of these systems is the existence of 
discrete analogue of equations which preserve integrability. This is remarkable since there 
are infinitly many possible ways of discretization which violate integrability and create 
chaos. We are interested in understanding criteria which distinguish integrable discrete 
systems from nonintegrable ones. In the case of 2nd order ordinary differential equations, 
Painleve property plays one of such roles. 

In this paper we consider a discrete analogue of the Painleve property and show 
explicitly how it works to characterize integrability. We study, in particular, iV-point 
Lotka-Volterra map. It is equivalent to the discrete time Toda lattice, hence is completely 
integrable. Such integrable systems are usually believed that their orbits are also deter- 
ministic and insensitive on their initial values. Namely starting from an arbitrary initial 
value an orbit is determined uniquely. 

The purpose of this paper is to show that above view of a completely integrable system 
is not always correct. By analysing the discrete Lotka-Volterra equation we will show that 
the orbit is not determined uniquely but there arise two ways of orbit to be chosen at every 
step of time evolution. This unexpected behaviour owes to the appearance of a square root 
function. Hence the situation is similar to the Painleve property for ordinary differential 
equations, in which an appearance of a moving singularity implies nonintegrability of the 
equation. We can show explicitly the branch point of the square root function moves 
depending on initial values when the equation is perturbed by a small parameter. If the 
parameter vanishes, the square root is resolved as expected, but two valuedness remains 
and causes uncertainty in its evolution. 

The equation we discuss in this paper is the following type of discrete time Lotka- 
Volterra (d-LV) equation : 

a£"(l-&£S) = 4(l-<^ +1 ), n = l,...,JV (1) 
Xq = x l N) x\ = x f N+1 (periodic condition), (2) 

where 5 is the minimum size of time step. 

In the limit 5 — > , ([!]) becomes the continuous time Lotka-Volterra equation : 

-^ x n(t) = x n (t)[x n _i(t) - x n+1 (t)}, (3) 
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which is more familier and known integrable. Conversely equation (|I|) is a natural dis- 
cretization of d||) in the sense that it is obtained from (||]) within the framework of Hirota 
bilinear formalism, which preserves integrability. 

The complete integrability of the equation (HD itself is guaranteed by the existence of 
sufficient number of conserved quantities. For our purpose it will be more useful to notice 
the connection of d-LV (P to the discrete time Toda lattice, 

I nit + 5)V n (t + 6) = I n+1 (t)V n (t), 

(4) 

I n (t) - I n (t + 6) = 6 2 [V n . l (t + 6)-V n (t)}, 
by the Miura transformation |l| 

V n (t) = x 2n x 2n+ ii 

(5) 

I n (t) = (l-fei^xi-^y. 

The discrete time Toda lattice is a completely integrable system which has been well 
studied. Many useful information can be translated from the Toda lattice to our system. 
At the same time all results we obtain for the d-LV in the following discussions must be 
true also for the discrete time Toda lattice. 

Recently integrable discretization has been further extended to ultra discretization. 
Namely associated to variety of discrete soliton equations there exist integrable equations 
whose dependent variables are also discrete. Starting from certain class of discrete soliton 
equations we can derive so called ultra discrete equations in a systematic way. In this 
way the ultra discrete analogue of d-LV was discovered in |J : 

d^ 1 -4 = max[0,4 +1 -l]-max[0,4t 1 i-l] , n = 1, 2, • • • , N. (6) 

We will also discuss this equation from the same point of view. 

We are going to examine in detail the behaviour of solutions of N point discrete Lotka- 
Volterra system in two separete papers. In the present paper we concentrate our attention 
to integrable case. Although solutions have been given explicitly in terms of r functions, 
we will show that an actual orbit is not deterministic. At every step of time evolution 
there are possibility to jump from one stable orbit into another. Two different orbits are 
characterized by the same constants of motion. This unexpected behaviour exists only in 
the case of discrete map. Therefore we will show how this uncertain behaviour ceases to 
occur in the continuous time evolution. We are also interested in whether both types of 
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solutions remain in the ultra discrete limit. In conclusion we will find that the orbit is 
determined uniquely. 

Our study of the discrete Lotka-Volterra system can be extended to nonintegrable 
case by introducing a small deformation parameter. We can see in detail behaviour of 
the singularities at every step of time evolution. In addition a Julia set will turn out to 
be a proper object to study global feature of singularities. We can see how the Julia set 
disappears from the complex plane of dependent variables in the integrable limit. We will 
discuss all these in our succeeding paper. 

2 3-point Lotka-Volterra system 

First we concentrate our argument to 3-point (N=3) case, but most of the results which 
we discuss below can be extended to arbitary N cases. Since we are interested in studying 
difference of variables before and after one step of time evolution, we set 5 = 1 and denote 

(x, y, z) := [x 1 , x 2 , x 3 ), (x ,y , z) := (x^~ , x 2 + , x 3 + ). 

Under this simplification of notation the 3-point d-LV equation is 



x'\ 


;i - 


= X[ 


;i 


-y), 


(7a) 


y'i 


l-x>) 


= y( 


i 


-z), 


(7b) 


z'i 


;i - v') 


= z( 


;i 


— x). 


(7c) 



Applying the scale transformation (x,y,z) — > (5x,5y,Sz), we can restore the minimum 
time lemgth 5 and recover the equation ([!]) for the 3-point case. Hence it is sufficient to 
consider only 5 = 1 case (7), unless the continuous limit is taken. 

2.1 Two types of map 

In order to see the time evolution of the system we solve (7) for (x', y', z') as functions of 
(x, y,z). For convenience we denote the right hand sides of (|7a|) , ( [75| ) and ([Tc|) by 

X = x(l-y), Y = y(l-z), Z = z(l-x). (8) 
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Then the successive substitutions of ( [7T5D and ( [7c] ) into ( [7a] ) yield 

, X X X . , 

x = l 17 = 5^ = v (9) 



i - y' y 

1 - x' 

i - y - x' 

^ „ 7, 7= (10) 



y - z - (i - z)x' 



Solving this for x' we obtain 



. i + x- y- z±v/(i-x-y-z) 2 - axyz . . 

X= 2(1 = *) • (U) 

From this elementary calculas it is apparent that the orbit of the map is not unique. 
The singularity of the map associated to the appearance of the square root function arises 
at every step of the map. Since there is no reason for a solution in (|TT] ) to choose one 
of two signs of the square root, the map is two valued, hence not unique, at every step 
of the time evolution. The uniqueness of the map is violated unless the argument in the 
square root vanishes by itself irrespect to the values of x, y, z. But our result flllf) shows 
that it vanishes only if x, y, z are in some relation. 

How this could be compatible with the complete integrability of the equation. In order 
to clarify this point, let us substitute (§) to ([TT|) and see what happens. First we notice 
that the square root is resolved, i.e., 



- X - Y - Z) 2 - AXYZ = \xyz - (1 - x)(l - y){\ - z)\, (12) 

hence the map is not singular. This property resembles the Painleve property of ordinary 
differential equations. To make this point clear we must compare above result with the 
case where X, Y, Z differ from (§). We will, however, postpone this argument to our 
succeeding paper. 

If we substitute this result into (p[) , we obtain 

X l-z + zx (13) 
1-2/- 

Hence we are left with the problem of non-uniqueness of the map. Both in ( ]TB| ) are 
candidates of solution to (7) with equal importance. At every step of evolution there is 
no reason to choose any particular one out of two solutions. If we repeat the map m steps 
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there will arise 2 m possible orbits which we can choose both with same weight. An orbit 
obtained by solving (7) seems to become non-deterministic. 

Now we ask how this could be compatible with complete integrability of d-LV. This 
will be understood if we solve (7) also for y' and z' in similar way as we did for x' case. 
Because d-LV equation has cyclic symmetry under the rotation of (x,y,z), the result is 



( i 

X- 




y 



y 



1 — Z + ZX 

1 — z + zx 
1 — x + xy 

1 — x + xy 



A- type 




5-type (14) 



V l-y + yz/ 

Hereafter we call these maps A-type and B-type respectively. 

We must notice that there are not eight but only two types of the map at every step 
of evolution. If the minus sign is chosen in (|TTP, the same sign must be chosen also for 
y' and z' to obtain a consistent answer. Do we still get 2 m different maps after m steps? 
This is not so because of the special character of the B-type solution. Namely it is a 
map of reflection. If we repeat the map six times a point goes back to the original place. 
Moreover two successive map of B after the A-type map is equivalent to the two successive 
map of B followed by A-type of map. We can write them symbolically as 



B k 



AB = B A. 



Similarly we can find many such relations. Therefore the number of possible orbits will 
be reduced to smaller, but not to one. 



2.2 Conserved quantities 

Now we consider relation between conserved quantities and the two types of the map. In 
the 3-point case, conserved quantities are easily found, 



C = X+Y+Z=x+y+z— xy — yz — zx 
D = XYZ = xyz(l-x)(l-y)(l-z) , 



(15) 
(16) 



which are obtained simply by summing or multiplying the left hand side of (7), respec- 
tively. There is no other quantity independent from these two and being constant of 
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motion. These two constants restrict orbits of the motion to some curves in the three 
dimensional real space (x, y, z). 

To see the orbits more explicitly it will be more convenient to introduce new constants 
r, s defined from C and D by 

rs = D 

r + s = 1-C. (17) 

Writing them explicitly we have 

r = xyz, s = (1 — x)(l — y){\ — z). (18) 

The constraint of r to r = a (const) forms a hyperboloid. Similarly s = b (const) forms 
another hyperboloid. Hence the orbits of the map are bound into intersections of these 
two hyperboloids. 

We now observe that if the initial condition (r, s) = (a, b) fixes orbits for a set of 
constants (C,D), (|i~7D shows that (s, r) = (a,b) also fixes orbits with the same values 
of (C, D). It amounts to exchange (x, y, z) with (1 — x, 1 — y, 1 — z). They form differ- 
ent intersections which are point reflection of each other with the center of reflection at 
(x,y,z) = (1/2,1/2,1/2). Therefore for a given set of constants (C,D), there are two 
sets of separate orbits, which we will call + and 0_. Figure [I] illustrates an example 
of these maps, which are projected on the x, y plane. In this figure initial values are 
chosen at (x,y,z) = (0.06,0.6,0.17) which draw 0+ with (r,s) = (0.00612,0.31208). On 
the other hand 0_ is drawn with (r, s) = (0.31208,0.00612) by starting from (x,y,z) = 
(0.94,0.4,0.83). 

We are now ready to see correspondence between the orbits drawn in the picture and 
our solutions in ([T4]) . By looking at the A-type solution, we find that (r, s) remain the same 
after the map. On the other hand a B-type solution exchanges (r, s) to (s, r). Needless to 
say both maps are allowed since the constants of motion C, D remain unchanged. Thus 
we obtain the following correspondence 

• A — type r — > r , s —>■ s (r,s are unchanged) 

• B — type r — ► s , s — ► r (r,s are exchanged) 

Starting from particular initial values of (x,y,z), say (0.06, 0.6, 0.17), and if there is 
no B-type solution, we will find only + orbit. 
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Figure 1: Example of orbit (time steps 300) 
2.3 Continuous time limit 

In section 2.1 we have seen d-LV eq. has two types of map. We are interested in their 
behaviour in the continuous time limit. The d-LV equation itself turns to the continuous 
time LV equation @ in the 5 — > limit. This fact seems rather paradoxical in the first 
sight, because the uniqueness of solutions must be true in the case of ordinary differential 
equations. Namely there is only one orbit starting from a particular initial value. We like 
to see how the correspondence between discrete system and continuous system holds in 
our case. We will try to understand this problem from two different view. 

First we look at the solutions (|HD in detail when 5 is small. For this purpose we 
restore 5 in fll3|), 

x i — x 1 1 x2 1 1 : A ~~ typ e 

1 — 0Z t + Z Z t X t /^g\ 

i - y* : B - type. 

Expanding the both sides of the A-type solution for small 6 we find 

x t + 5-- + o(5 2 ) = xHl - <V + 6z* + o(S 2 )) 
at 

Hence in the first order of 5 we obtain the continuous time LV equation (|3|). The same 
arguments are true for y* and z*. This means that the A-type solutions satisfy the LV 
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Figure 2: Change of orbits under variance of parameter 5 

equation (3) in the continuous time limit. As for the B-type solution all points are mapped 
to infinity when 5 turns to zero. In other words the B-type solution has a meaning only 
if S is finite. In this sense the solution becomes unique in the continuous time limit. 

The second view to clarify the behaviour of solutions is to see the dependence on 5 of 
the cross section of the hyperboloids specified by the constants of motion. We have seen 
that there are two different kinds of cross sections, which we called 0±. Under a variance 
of 5 we will find that the orbits 0_ move far away to infinity as 5 becomes small, whereas 
the other orbits + remain in finite region. 

Figure shows such an example of change of orbits for 5 = 1.0, 0.8, 0.6 respectively. 
Initial values are fixed at (x,y,z) = (0.06, 0.6, 0.17). 

3 N-point Lotka-Volterra system 

We can generalize most of the results for 3-point discrete Lotka-Volterra system which we 
found in the previous section to arbitrary iV point case. 
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3.1 Two types of map 

The equation which we concern in this section is (|l|) with arbitrary N. Let us simplify, 
for a while, the notation by putting 

x' n = 6xi +5 , x n = dxi (20) 

and the equation of motion 

x' n {l-x' n _ l )=x n {l-x n+1 ), n = 1,2,..- ,N. (21) 
We also denote the right hand side of (|2~lD by 

X n :=x n {l-x n+1 ), n = 1, 2, • • • , JV, (22) 
and solve (|2l| ) for x[ in terms of (Xi, X 2 , • - - , Xn)- Repeating the same procedure as in 



the case of 3-point d-LV we find 



x\ = % (23 ) 

j AM 



X 2 

1 - 



1 - x[ 



= Xl c *-i- a *-fi (24) 

Cn — Qjv^i 

where a^, satisfy the recurrence relations 

a k+1 = a k -a k -iX k +i, (a , eti) = (0, 1), (25) 
Cfc+i = c fc -c fc _iX fc+1 , (c ,ci) = (1, 1), (26) 

or solving them explicitly 

a k = a(3,fc), c fc = 0( 2) fc) (27) 

where 

„, r _VrW) / -/ [( fc -^ + 2 )/ 2 ] (OH) 
°(i,fe) — 2^°» ' W-f + £ + 2)/2] fc<z 
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and 



a 



1, c x 



m ._ 



(i,k) 



-iy J2 x h x i2--- x i» i*o- 



J1,J2,--- ,Jl 



The summation ^2^% ■■■ j l mus t De taken over j a 's satisfying 



J1J2, ■■• JiE 



i < fc 



(M + l, ■■• ,N, 1,2,- ■■ ,k), i>k 

\ja~Jb\>2, 

and [k] is the maximum integer which is equal or smaller than k. 
From (]24]) becomes 

, b ± a/& 2 — 4a A rCiv_iX 1 



2a 



6 := + ajv-i^i- 



A' 



If we define d by 



d := a/& 2 — 4a iV CAr_ 1 X 1 
and use the relation (|69|) in the Appendix for i — 1, we obtain 

6 = c^v — a A r — 1-^1 4" 2ajv_iXx 

[7V/2] 

= Ci + 2a A r_ 1 X L , 

z=o 



where 



(1,A0 

Cl ■= C l ' = (-1)' ^ X h X h ' ' ' X or 

31,32,— ,31 



Moreover, using the result of Appendix (73), 



{cncln-i — <^nCn-i)Xi — XiX 2 ■ ■ ■ X^ 



we find the expression for d 



{[N/2] \ 2 

^E^l -4A 1 X 2 ---A 



AT] 



«=0 



which is symmetric under the permutation of X n 's. 

We have arrived at the expression which looks complicated. The map could be 
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valued only if d vanishes and could be regular only if the square root is resolved. How 
the equation could be completely integrable. 

A miracle happens if we substitute (p2|) into ( ]34l) . The square root is resolved and we get 
(see Appendix) 



d 



x\x<i • • -Xjsi — (1 — x\){\ — x 2 ) •••(! — xn) 



(35) 



At the same time we also observe that the two-valuedness of the solutions is inevitable. 
Therefore the orbits will split to two at every step of the map. The situation is just 
like it happend in the 3-point case, but the expression of solutions looks much more 
complicated. This is, however, not a problem, because we can find one solution out of two 
easily. Namely if we suppose a set of solutions for all n and substitute 

into the equation ([H]), we will immediately see that it is satisfied. Once we know one of 



solutions, say X-, the other one is readily found to be b/a^ — X- from ([H|). In summary 
we have obtained 

b 

l + x 2 

x\ = { aN . (36) 

1 — x 2 

We obtain general solutions from ( |36l) by symmetric permutation of arguments: 

2)X n ±d 

*n = 1 1, " ] " I 37 ) 

/a (n+2,n-l) 

fl(n+l,n-l) + a >{n+2,n-2)X n 

1 + X n +i 



0(n+2,n-l) 



(3? 

1 — x n+l 



0-{n+l,n-2) . 

x n A - type 

0( n +2,n-l) 



(39) 



1 - x n+1 B - type 

To obtain the last expression, ([74]) is used. 

The result we have obtained is quite similar with the case of three point d-LV equation. 
There are two types of distinct map at every step of time evolution. The B-type is a map 
such that B 2 is a cyclic permutation of points. When is even B N = 1, and when A^ is 
odd B 2N = 1 holds. We also find that AB 2 = B 2 A is satisfied. 
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3.2 Conserved quantities 

We like to show, in the following, that the CVs which were defined in (|32|) are constants of 
motion. To prove that we substitute (p^ ) to (|32|) and expand factors and then recombine 
terms of the same degree of powers. We obtain 

2j ^ji^j'a ' ' — — X j2+ i) •••(!— Xj l+ l) 

ji,h,— >ii 
1 (i.*0 

= X 3l X 3l' ' ' X 3l ^JO! +l' r ia 2 +l ^jar+l 

r=0 jl,h,— ,3l 0,1,0,2,-, a r 

r=0 ji,32, ai,a2,— ,a T 

To obtain the last expression we have used the fact that, owing to the cyclic symmetry 
of the summations, there are always the term Xj a -\Xj a -1 ■ ■ ■ Xj ar -i in the summation if 
there is Xj +\Xj +1 ■ ■ ■ Xj ar+ \. Then we can recombine terms to get 

^ ] x h x h ' ' ' x ji(^- ~ x ji-^)(^ ~~ x 32—l) " ' ' (1 — x ji-l)- 

ji,h,— J 1 

Now by usig the equation of motion (f2~l~D, this is nothing but Ylji j 2 ... j, -^-31-^-32 ' ' ' -^ji 
evaluated at time i — 1. Therefore we conclude that C\ with Z = 1,2, ••• , [iV/2] are 
constants of motion. 

Now we consider relations of the two time-maps and the conserved quantities. To this 
end we first notice that X\X 2 ■ ■ ■ X^ is also conserved, as it is apparent from the equation 
of motion. Since C/'s are constants, d of (pi) is a l so a constant. If we define 



r := X\X 2 - ■ ■ x N , s := (1 — Xi)(l — x 2 ) ■ ■ ■ (1 — x N ), (40) 

we have 

rs = X\X 2 ■ ■ ■ Xjsr, \r — s\ — d. 
Hence we find that rs and |r — s\ are both constants of motion. 

In the 3-point d-LV the B-type maps are characterized by the exchange of r and s, while 
they remain unchanged by the A-type maps. This is true also for arbitrary iV point d-LV. 
For example we see, for the A-type map of (|39D , 

1 1 1 a (2,-l) a (3,0) a(N+l,N-2) 

X \ X 1 ' ' ' X N — X\X 2 ' ' ' X • • ■ 

a (3,0) a (4,l) a(N+2,N-l) 

= xix 2 ---x N 
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Since rs is an apparent constant this result also means s does not change. Therefore we 
conclude 



• A-type r — > r , s — > s (r,s are unchanged) 

• B-type r — > s , s — > r (r, s are exchanged) . 

3.3 Continuous time limit 

From the construction of d-LV equation there is no problem to obtain (H) as the continuous 
time limit of ([!]). As it was the case of 3-point d-LV, this does not mean that all solutions 
of d-LV remain being solutions of ([I]) in the limit 5^0. 

For the A-type map of (|39D , collecting the terms of lowest order in 5 we obtain 

n + dt + ° {d) ~ X - 1 - 5{xi +2 + • • • + + o(^) (41 j 



or 



dx 



dt 



t 

n (ft I . . . I rpt I ryX 

~ • h n\- b n+2~ ' n-\-N—2 ' n-\-N—\j 



X n( X n+l + X n+2 + ' ' ' + X n+N-2) 



X n( X n+l X n-l)- (^2) 



Hence we find that the A-type map of ( ]3"5| ) reproduces (3), which means that every solution 
consisting of a sequence of A-type map is a solution of (^) in the continuous time limit. 
As for the B-type map of (PD|) we have the expression 



_ 1 _ t 



after restoring 5. All the solutions of this type tend to infinity in the limit 5^0. 



4 Ultra discrete limit 

The ultra discrete Lotka-Volterra equation is another integrable deformation of the LV 
equation. Since it is derived starting from the d-LV equation which include both A-type 
map and B-type map, we are interested in how the solutions to d-LV equation behave in 
this limit. 
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4.1 3-point case 



The ultra discrete version of LV equation was shown to exist in || . 3-point ud-LV equation 
under periodic condition is: 



d[ — d\ = max[0, d 2 — 1] — max[0, d' 3 — 1], 
d>2 — d% = max[0, d^ — 1] — max[0, d\ — 1], 
d's — c?3 = max[0, g?i — 1] — max[0, d' 2 — 1], 



(43) 



where max [a, 6] means to choose a if a > b and 6 if 6 > a. This equation is derived from 
such transformation of variables: 



5 = exp[-l/e], 



X r 



-exp(c4/e) 



(44) 



and ultra discrete limitQ e — > +0. 

First let us attempt, like the d-LV case, to solve ud-LV equation flUf) for (g^, e^, in 
terms of (g?i, g? 2 , d 3 ). Successive substitutions yield 



d\ + max[0, c/ 2 — 1] — max[0, d' 3 — 1] 
d\ + max[0, d 2 — 1] 

, g?3 + max[0, g?i — 1] — max[0, d' 2 — 1] 

d\ + max[0, d 2 — 1] 

, g?3 + max[0, d\ — 1] 



— max 



max 



max 



, 6?2 + max[0, c/3 — 1] — max[0, d[ — 1] — 1 



-1 



which is to be solved for d[. The appearance of d[ in max statement in the right hand 
side prevents us to proceed this calculas further. Thus we must seek other method to find 
solutions. 

For the meantime we suppose (di, d 2 , ofa) have a specific set of integers. Then we can seek 
(d[, d' 2 , d' 3 ) following to the prossesses below: 

1. Suppose d[ is an integer a. 



When we derive 



we also use the property 

lim elog(l + e A7e ) =max[0,X]. 

e— >+0 V 
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2. From the value a and the second relation of (E3p we determine d 2 . 



3. Similaly we can determine d 3 from above d 2 and the third relation of ([43]). 

4. Check if these values (d[, d' 2 , d' 3 ) satisfy the first relation of (f43"|) . If so, (d[, d' 2 , d 3 ) is 
a correct set of time evolution. 



5. Return to 1. and try another integer. 

One may think there are more than one correct set of time evolution. But we find 
numerically only one correct set of (d[, d 2 , d 3 ) starting from a given set of (di, d 2 , d^). We 
have checked this for every integer a between ±1000. In this way we are convinced by 
ourselves that (d[, d 2 , d 3 ) is uniquely determined. 
Figure 13] shows some examples of such ud-LV system. 



t = 
t = 1 
t = 2 
t = 3 



(2 
(1 
(1 
(2 



1) 
2) 
1) 
1) 




1 

2 
3 
4 
5 
6 



2) 

1) 
0) 
0) 
1) 

2) 
2) 



t = 





(8 





4) 


t = 


1 


(0 


3 


9) 


t = 


2 


(2 


10 


0) 


t = 


3 


(11 





1) 


t = 


4 


(1 





11) 


t = 


5 


(0 


10 


2) 


t = 


33 


(8 





4) 



Figure 3: Example of 3-Point ultra discrete Lotka-Volterra system 



Now let us see ultra discrete limit of solutions. First we will consider ultra discrete 
limit of A-type. We again start from the A-type solution of (|19|) . By the transformation 
(U), it turns to 

M_ ldi l + e^- 1 > + ei^ +d 3-2) 
6C 1 CE i + e |(*-i) + e i(d3+d!-2) ' 

hence 

d[ - d x = eln (l + e^" 1 ) + e §(<fe+<*3-2)^ _ eln ^ + e |(*-i) + (4g) 
The ultra discrete limit e — >■ +0 yields 

d[ = d x + max[0, d 2 - 1, d 2 + d 3 - 2] - max[0, <i 3 - 1, d 3 + ^ - 2]. (47) 
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Since there is no d[ in max statement in the right hand side, we can get directly information 
for d[ from preceeding ones. The general formula for other d n can be found from ( fl7j) as 

d' n = d n + max[0, d n+1 - 1, d n+1 + d n _ x - 2] - max[0, d n _ x - 1, d ri _i + d n - 2] (48) 



From fl48|) we find immediately the following results: 

• If all of d n are less than 1, all d n remain constant. 

• If all of d n are greater than 1, d' n — d n _\. 

We will have more complicated pattern of time evolution if there are <i n 's, some of which 
are less than 1 and some others are greater than 1. 

As for the B-type map the solution in (^) becomes, after the transformation ([H]). 

l + e M-i) + e *<*-i) = o. (49) 

By taking logarithm of both sides and letting e go to +0 we get 

lim eln (l + e^i" 1 ) + e~S d ^\ = maxfO, d[ -l,d 2 - 1] (50) 

= ex(-oo). (51) 

This is consistent for none of <i„'s. Therefore we conclude that the B-type map does not 
give any meaningful result in the ultra discrete limit. 

4.2 Generalization to N-point case 

iV-point ultra discrete LV under periodic condition has been given byH 

d 1 ^ 1 -d t n = max[0,d t n+1 -l]-max[0,d t ^_ 1 l -l], n = 1, 2, • • • , N. (52) 



This can be derived from (|l|) using the rule of transformations (j44[) . 
By the same reason as in the 3-point ud-LV case we cannot derive simple ud-LV time- 
map directly from (|52|), since the unknown variables are under the max statement. We 
can, however, derive it from the A- type map (|39|). According to the prescription of ultra 
discretization (fH|) and taking the ultra discrete limit e — > +0, we have 

d' n = d n + lim (elna (n+1)n _ 2 ) - elna (n+2 ,n-i)) 



17 



The ultra discrete transformation of Xj is given by 

Xj = 5xj(l — 5xj + i) —exp((dj — l)/e) — exp((dj + dj + \ — 2)/e) 



Since, in a( n+ i in _ 2 ) and a( n +2,n-i), -X/s appear in the form (— VfX^X^ ■ ■ • Xj v 

[([fc-i|+2)/2] (i,k) 



Hhk) 



ud 



Yl exp +d h + --- + d, h -l)/e) 

1=0 ji,h,— ,h 



x {1 + exp((d il+ i - l)/e) + exp((rf i2+ i - l)/e) + • • • + exp((d ji+1 - l)/e) 

+ exp((dj 1+1 + d j2+1 - 2)/e) H h exp((d il+ i + d, 2+ i H h - Z)/e)} . 

Let us denote by (a, 6) © (c, d) the set (a + c, a + d, b + c, b + d). Then we have 



lime(ln a^)) = max 



[flfc-i|+2)/2] 

(J (J © (m ix © m j2 © • • • © m A ) 



where 

TOj := (dj — 1, dj + cfj- + i — 2). 
Finally we obtain the solution of ud-LV as 



d' n = d n + max 



— max 



[(AT-l)/2] (n+l,n-2) 

|J |J © (m^ © m i2 © • • • © m y J 

i=o h,h,— ,5i 
'[(JV-l)/2] (n+2,n-l) 

U U ©K®^©--©^) 

«=0 h,h,— ,3i 



This map will give the same map derived from (52). It, however, comes only from the 
A-type map. As for the B-type map we will not get a meaningful answer by the same 
reason we discussed in the case of 3-point ud-LV equation. We could guess this result 
because the ultra discrete limit also implies the continuous time limit as we could see 
from (PI). 



5 Concluding remarks 

In this paper we have discussed properties of discrete Lotka-Volterra equation. A special 
attention has been focussed on the fact that the map is not single valued at every step of 
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the time evolution. This singular behaviour disappears both in the continuous time and 
ultra discrete limits, and uniqueness of the map is fulfilled. 

The violation of the uniqueness of the d-LV system does not contradict with complete 
integrability, in the sense that it has sufficient number of constants of motion. A singular 
map, which we called B-type, does not create new orbits during the evolution but simply 
causes jumps of the system between two possible orbits. 

We have explained in detail how these two maps arise in the d-LV system. Usually such 
a solution as B-type was not taken into account to study a deformation of the continuous 
time Lotka-Volterra equation ||. In principle, however, there is no a priopi reason to 
ignore the B-type solutions. 

General solutions are random sequence of A-type and B-type maps. In other words 
two types of solutions can be chosen arbitrarily at every step of the map, hence we can 
not write their orbits in a compact form. If only A-type of map is selected every time 
many solutions are known explicitly by means of the r functions of the theory of Toda 
lattice. We will present here a particular one of such solution which corresponds to the 
cnoidal wave solution of the continuous time Toda lattice, i.e., 

E 



V n (t) = (2Ku) 



(53) 



where dn is a Jacobi elliptic function, K and E are called 1st and 2nd complete elliptic 
integrals, respectively, and v and A are related by the spectral relation 

In terms of r functions the function V n of the Toda lattice is expressed as 

r'+V 

W) = J ^r ( 55 ) 

n n 

hence, by means of the Miura transformation (|5|), the function of d-LV is given by 

< = r( ^r r H i)/2 - (56) 

T (n-l)/2 r n/2 

Corresponding to the cnoidal wave (|53| ) r* is given by 

Ti = 0t>(yt±^ + t) (57) 
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where $0 is the Jacobi theta function and is an arbitrary constant phase. Substitution 
of this into (|56|) yields 



# ( v (t + l) ± + <i>)0 o (rt ± *g + <!>) 



- A (59) 

where 

f=2ir(w±£ + *) , i = k( v *±) = (60 ) 



and is the modulas which characterizes the elliptic functions. We have used 
1 - fc sn f sn 77 = , f = 2Kx, r] = 2Ky. 

#oOw(y) 

The Jacobi sn function has a symmetry 

sn(£ + 2K) = -sn£. 

Comparing (|59D with @) the periodicity condision on x l n requires for A to satisfy 

N 

A — — , I — 1, 2, 3, • • - . 

In the case of three point d-LV, A = 3, 3/2, 1, 3/4, • • • . 

If only B-type of map is allowed to exist, we have general solution in the following 
form: 

{/(n + t) iit odd 

1 (62) 
- — f(n + t) if t even 



where f(m) is an arbitrary function of period N, i.e., f(m + N) = f(m). 

It will be interesting to understand what is the corresponding behaviour of the B-type 
map in the discrete time Toda lattice. The answer is rather simple. It corresponds to 
exchange the role of V and J, as we can see directly in the Miura transformation (|5]). This 
property of the Toda lattice is transformed into two distinct orbits, which we called 0±, 
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of the d-LV system. We note that this symmetry is different from what is called duality 
under the V — I exchange in the case of continuous time ||, since our symmetry arises 
only in the discrete time evolution. 

Before closing this paper we like to point out that the singular behaviour of the map, 
which we discussed, is not special in d-LV equation. As we will explain in the other 
paper, this nonuniqueness of the map turns out to be a source of Julia set which should 
be generated when the equation of motion is deformed from the integrable case. In order 
to make this point more clear we recall how the integrable logistic map is changed to the 
non-integrable one. 

In a series of our previous papers || we have studied a generalization of logistic map 
which interpolates between the Mobius map and the logistic map. The former is integrable 
whereas the latter is well known non-integrable map which generates chaos. To be specific 
we consider 

Z t+1 =fMZ t — . r— . 63 

1+^(1-7)^ 

7 = and 7 = 1 are the Mobius and the logistic map respectively. The Mobius map can 
be integrated explicitly as 



a + jj,'- 



« = -^j^. (64) 



1-f 

If 7 7^ 0, the map ( p5| ) is not integrable. Depending on its initial values as well as on 
the values of parameters /1 and 7, orbits converge to some attractors, otherwise behave 
chaotic. In order to see if the map is integrable or not, it is more appropriate to study 
Julia set of the map. The Julia set is defined on the complex plane of z t . It is a closure 
of repulsive fixed points, which can be obtained by the inverse map starting from an 
arbitrary repulsive fixed point. Since it is an invariant set of the map it does not depend 
on which repulsive point is chosen at first. There exists a Julia set on the complex plane if 
the map is not integrable. Therefore we are interested in how the Julia set will disappear 
in the integrable limit 7^0. 
The inverse map of is readily solved 



_ 1 - (1 - j)zt +1 ± y/(l - (1 - 7)^+i) 2 - 47*t+i//x 
z t - - . (65) 
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This map is not unique and the Julia set has a fractal structure. In the integrable limit 
7 — > we obtain 

z t = { Ml - Zt + l] (66) 

hence the square root singularity is resolved but the map is not unique. The 'B-type' map 
does not create new orbits. It maps every point to the origin. If we produce the Julia 
set by using computer we will not observe this point. But it is not right to forget this 
point. When the equation is deformed from the integrable one two maps of ( |65D are equal 
footing and become the source of chaos. The situation is exactly the same with what we 
learned in this and the following papers. 
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Appendix 



Let us start from the definition (|2q): 

^ [(.V - /' + /,- + 2)/2] < i 



j '"max 

1=0 

This satisfies the recurrence formulae 



j _ J [(fc-z + 2)/2] k > i 

1=0 



a (i-l,k+X) — a (i-l,fc) — a (i-l,fc-l)^fc+l ? (67) 

= a>(i,k+i) — a(i+i,k+i)Xi-i, (68) 



when k — % ^ — 1, —2 mod N, and 



a (M-l) — — a (i+2,i-2)Xi, 

= 0(i,i-2) — a>(i+l t i-3)Xi-h 5 (69) 



when k — i = — 1 mod N. 
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Using pi) 



Similarly using 



A 



i k 



a (i,k) — (1 — ^i) a (i+l,fc) 

a {i+l,k) — a {i+2,k)Xi — (1 — Xj)a(i + x,k) 

XiAi+i^ — • • • 



B, 



a (i,k) ~ %k+l a (i,k-l) 

a (i,fc-l) — a {i,k~2)Xk — Xk+ld(i,k-l) 

-Bj,fc-l(l — JEfc+l) = ■ • • 

(1 -rr<) • • ■ (1 - x k )(l -x k+l ). 



Let us prove (pq). Using (p9|) 



[JV/2] 



i=0 



C« — — Q(i,i-2) — Q(i+l,i-3)^i-l 

— ^i,i-3(l — x i-l) + ^-i,i-Z x i-l 
N N 



3=1 



3|) can be shown by using ( j6^) as 



a(i-l,fc+l) a (i,fc) — a {i-l,k)^(i,k+l) — [a(i-l,k) a (i,k-i) ~ a(i-l,k-l) a (i,k)) 



-Xi_iXi ■ ■ ■ X k+ i 



Finally we prove fl39|): 



{Xn+l ~ l)0(n+2,n-l) + a (n+l,n-l) + 0( n+ 2,n-2)^n 
^4n+l,n-l — A n+ i, n -2X n + fl(n+l,n— 2)^n 
a (n+l,n— 2) a 'n 



where (1701) is used. 
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